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Abstract 

We study the Hamiltonian formulation for a parametrized scalar field in 
a regular bounded spatial region subject to Dirichlet, Neumann and Robin 
boundary conditions. We generalize the work carried out by a number of au¬ 
thors on parametrized field systems to the interesting case where spatial bound¬ 
aries are present. The configuration space of our models contains both smooth 
scalar helds defined on the spatial manifold and spacelike embeddings from the 
spatial manifold to a target spacetime endowed with a fixed Lorentzian back¬ 
ground metric. We pay particular attention to the geometry of the infinite 
dimensional manifold of embeddings and the description of the relevant geo¬ 
metric objects: the symplectic form on the primary constraint submanifold and 
the Hamiltonian vector fields defined on it. 


Key Words: Parametrized field theories; bounded domains; Hamiltonian formula¬ 
tion. 


1 Introduction 


Parametrized theories are interesting examples of diff-invariant field systems. They 
were introduced by Dirac [1] circa 1950 in order to explore the idea of forgoing the 
use of hxed spacetime foliations with flat spatial slices, and consider a more gen¬ 
eral approach to the study of relativistic held theories. They were further analyzed 
by other authors, in particular Kuchaf, Isham, Hajfcek and Torre [2-5], in the case 
where the spatial manifolds were taken to be closed (i.e. compact without bound¬ 
ary). Parametrized scalar helds have been recently used to explore the resolution 
of a number of technical issues related to the treatment of diffeomorphisms in loop 
quantum gravity [6,7] and in the search for boundary observables [8]. Moreover, some 
interesting gravitational models, such as Einstein-Rosen waves in vacuum or coupled 
to massless scalar helds, are known to be dynamically equivalent to a parametrized 
massless scalar held with cylindrical symmetry [9,10] (see also [11] and references 
therein). 

The main idea in the parametrized framework is to adjoin space diheomorphisms 
as dynamical variables to the usual set of helds. This is done in such a way that 
the solutions to the standard held equations can be naturally mapped to the one of 
the parametrized model and viceversa. However it is important to point out that 
the interpretation of parametrized models dihers in some important aspects from the 
standard ones from which they are derived. One reason for this is the presence of 
gauge symmetries associated with the reparametrization invariance introduced by the 
diheomorphisms. 

The standard Hamiltonian formalism of mechanics is essentially geometric: dy¬ 
namics takes place in the cotangent bundle associated with a space of physical conhgu- 
rations, a diherentiable manifold itself, which is inhnite-dimensional in the case of held 
theories. It is then necessary to be aware of the difficulties that this dimensionality 
issue may present and choose the right mathematical tools to deal with them. As we 
will show in the paper, this is not difficult once the proper formalism has been identi- 
hed. A large part of the physical literature on parametrized held theories —certainly 
full of useful insights at the geometric level [12,13]— has some technical/practical 
disadvantages: the hrst one is the local character of the standard approach (use of 
particular coordinate systems and a certain forgetfulness about global issues), the 
second is the use of distributional objects at a basic level (i.e. the dehnition of the 
fundamental Poisson brackets). Although distributions can, of course, be dealt with 
and prohtably used in a completely rigorous way, they may pose some problems in 
practice, in particular for the type of systems that we want to study where man¬ 
ifolds of diherent dimensions (spatial manifolds and their boundaries) coexist in a 
non-trivial way. In this paper we will hence use coordinate-free, global, geometric 
methods. As we hope to show, once the proper set up is identihed, computations are 
straightforward. 

A very important feature of parametrized models is their singular character: their 
Lagrangians do not lead to a one-to-one hbre derivative (used to dehne the Legendre 
transformation). This means that it is not straightforward to obtain a consistent 
Hamiltonian description for them. The most popular approach to this problem relies 
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on the Dirac theory of constraints [14], A more geometric way to do it uses the 
so called GNH method developed by Gotay, Nester and Hinds [15-17]. Among the 
strong points of the latter are its global character and the possibility of naturally 
taking into account the necessary functional analytic issues relevant for held theories. 
In the particular case of parametrized models this method is very simple to use and, 
hence, will be employed here. 

The introduction of boundaries is interesting from a physical point of view. For 
obvious reasons, they are relevant in condensed matter systems but also in other 
contexts, specially in general relativity, where boundaries play a very signihcant role 
(spatial inhnity, black hole horizons, holography, etc.). A covariant Hamiltonian 
approach to general relativity coupled to different types of matter helds —in particu¬ 
lar, scalar helds— in bounded spatial regions is discussed at length in [18,19]. In the 
present paper we rely instead on the standard Hamiltonian description in phase space 
to study parametrized scalar helds. This is the simplest parametrized held system 
although some more complicated models can be considered (see, for instance, [20]). 
Our results can be compared with those of [18,19]— and checked to be compatible—by 
considering hxed foliations and non-dynamical metrics respectively. Notice, however, 
that our description uses the canonical phase space and, hence, the conditions in¬ 
volving momenta are not obvious in the covariant setting. An interesting application 
of the results that we give in the paper is the extension of the polymeric quantiza¬ 
tion carried out when no boundaries are present to these richer models with a more 
complicated phase space description [6,7]. 

To the best of our knowledge, reference [8] is the hrst paper which considers the 
Hamiltonian formalism for parametrized scalar helds subject to different kinds of 
boundary conditions (in particular of the Robin type). In that work the Dirichlet 
and Neumann cases are dealt with by resorting to standard Hamiltonian methods. 
On the other hand the Robin case is studied by using covariant Hamiltonian tech¬ 
niques because the boundary term introduced in the relevant action gives rise to some 
difficulties when going to phase space. As we show in the present paper the geometric 
and global methods that we use here can easily handle the Robin case while providing 
a very natural and simple interpretation for the results corresponding to Dirichlet and 
Neumann boundary conditions. We would like to mention here that our point of view 
is useful to understand in a clear way the bifurcation phenomenon in the constraint 
algorithm characteristic of parametrized gauge systems discussed in [20]. 

The structure of the paper is the following. After this introduction we discuss 
in section 2 the action principle for the parametrized scalar held with the boundary 
conditions considered in the paper: Dirichlet and Robin (of which Neumann is a 
subcase). Section 3 is devoted to the 3 + 1 decomposition of the objects that we use 
throughout the paper. The most relevant results related to the conhguration space 
and its tangent and cotangent bundles, where embeddings play a central role, are 
discussed in section 4. Section 5 is devoted to the Hamiltonian formulation. We end 
the paper in section 6 with a discussion of its main results and an appendix where 
we collect the variations of a number of geometric objects dehned on the manifold of 
embeddings. 
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2 The parametrized scalar field in bounded spatial 
domains 

Let us consider a four dimensional Lorentzian manifold (M,g) diffeomorphic to a 
product manifold [^ 1 ,^ 2 ] LI, where S is a smooth, orientable, compact, 3-manifold 
with smooth boundary (possibly empty). In this situation dM can be written as 
the union diM u dj^M u d 2 M where is diffeomorphic to [^ 1 ,^ 2 ] dH and diM 

{i = 1,2) are diffeomorphic to E. We will assume that diM are spacelike and that the 
Lorentzian metric g, with (e, +, +, +) signature, induces a Lorentzian metric gg ■= j^g 
on d^M, where jg denotes the pullback under the inclusion map jg : d^M ^ M. We 
have introduced the parameter e = -1 to allow for a straightforward extension of our 
results to the Euclidean (Riemannian) case. We will restrict ourselves to the case 
where {M,g) is time oriented. 

The starting point to arrive at the action for the parametrized scalar held that we 
study in the paper is 

f {^9~\dip,dip)-m^ip‘^)vo\g-^ f RV^volgg, (2.1) 

Z JM .Jd^M 

where (p : M ^ M is a real scalar held on M, R : M is a hxed smooth function 

and we are using the metric volumes on M and d^M. This action is the obvious 
generalization of the one used in [21] to discuss the Hamiltonian formulation for the 
scalar held with Robin boundary conditions (it can be essentially found in [22, page 
227]). Our results can be trivially extended to the case where a potential term R((p) 
is included as in [3]. 

The diherent types of boundary conditions that we consider in the paper are: 

1. If E is closed then we have the simple example of a massive (Klein-Gordon) 
scalar held without boundary. 

2. If E has a non-empty boundary and we demand (pojg = 0 we describe a massive 
scalar held subject to homogeneous Dirichlet boundary conditions. 

In these hrst two cases B plays no role and can be taken to be zero. 

3. If E has a non-empty boundary and we take R = 0 then we have a massive 
scalar held with homogeneous Neumann boundary conditions. 

4. If E has a non-empty boundary and we take R 0 then we have a massive 
scalar held subject to homogeneous Robin boundary conditions. 

In the last two cases no conditions on cp o jg are imposed a priori. Notice, also, that 
the Robin conditions include the Neumann ones as a particular example. 

Several comments are in order now. First, just mention that the non-homogeneous 
case is straightforward once the homogeneous one is understood (see [21]), so we will 
only consider the latter. A second comment refers to the way the variational principle 
works. We get the held equations by demanding that the action is stationary under 
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admissible variations of the fields with their values kept hxed at diM and d 2 M. The 
variations of the action are directional derivatives which may have both interior and 
boundary contributions. The presence of the latter does by no means imply that the 
action ceases to be differentiable. In fact, boundary contributions to the variations 
are actually the origin of natural boundary conditions such as the Neumann or Robin 
ones. 

Dirichlet and Robin conditions are implemented at the action level in different 
ways. In the hrst case one has to restrict oneself to held conhgurations (and, hence, 
also variations) that vanish at the boundary d^M, whereas in the second the boundary 
conditions themselves appear as a consequence of the requirement that the action is 
stationary under arbitrary variations. The Hamiltonian description of these models 
based on the use of the GNH method can be found in [21]. 

In order to parametrize the action (2.1) we introduce diffeomorphisms Z : / x E ^ 
M such that, for every t € I the image Z(Si) of the slice := {t} x S is 

spacelike and TZ.dt is timelike and future pointing. These diffeomorphisms are taken 
as dynamical variables in addition to the scalar held. We consider now the action 

f f {B‘^oZ)'il;^vo\z*gg, (2.2) 

Z JlxJ] Z JIxdT. 

where Z*g denotes the pullback of to / x S and '^:/xS-^Risa new scalar held. 

Regardless of the considered boundary conditions, the evolution of the scalar held 
^ is given by the Klein-Gordon equation 

{euz*g+m^)'ilj = d, (2.3) 

where the d’Alembert operator is dehned in terms of the Levi-Givita connection as¬ 
sociated with Z*g. There are no extra conditions on the helds if S is closed. In the 
case of imposing Dirichlet boundary conditions we must require that ^ vanishes at 
I X and in the Robin (Neumann, R = 0) case the variations of the actions provide 
the equations at the boundary 

(£TZ“hV(^) + (R2 o Z)V') =0 at J X as 

where TZ~^.V is the push-forward by Z~^ of the unit, spacelike, outer normal V to 
d’sM. No independent held equations are obtained by varying in Z as in the case 
when no boundaries are present, however it should be noted that the diheomorphism 
Z must take / x aS to a^M. 

3 Some remarks on the 3+1 decomposition 

We give here some results that we use to decompose geometric objects in / x S and 
get the Lagrangian for our system. To this end we introduce a unique decomposition 
of any tensor held on J x S with the help of the projector H := Id - df ® a^, which 
is a (l,l)-tensor held. For instance, vector helds Y € X(/ x S) can be written as 
Y = df(y)at+ n(y) and one forms a e Q^(I x S) as a = Q;(at)df + n(Q!). We will also 
use the notations H^ : X(/ x S) X(/ x S) and H* : f2^(J x S) ^ Q^(I x S) to refer 
to these maps. 
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Metrics 

Any metric g on / x S (in particnlar the pnllback g = Z*g that we nse throughont 
the paper) can be decomposed with the help of the identity Id = dt ( 2 ) 9* + 11: 


g = g{dt,dt)dt ®dt + dt 0 + g(<9t, If*-) 0 dt + g(n*-, If*-). (3.1) 

This is usnally written in the form 

g = (eiV^ + /3))df iS)dt + dt^P + ^^dt + j, (3.2) 

where the so called shift ft e {IxJ2) satishes ft{dt) = 0 , 7 is a symmetric ( 2 , 0 )-tensor 
satisfying j(-,dt) = 0 and 7 “^ (a symmetric (0,2)-tensor) is the nniqne psendo-inverse 
of 7 snch that 7 “^ ■ 7 = If and 7 “^(-,dt) = 0. Finally the nowhere vanishing fnnction 
N 6 C°°{I X S) is the lapse. The inverse metric g“^ can be conveniently written as 

g-' = •)) ® •))+r'. (3.3) 


Diffeomorphisms and embeddings 

Let Z:/xS^Mbea diffeomorphism snch that c M is spacelike for each 

f 6 /. We can dehne a decomposition of Z and its tangent map TZ adapted to the 
product structure of / x S by using the projector tensor held 11 introduced above. 
The tangent map TZ : T{I xT) ^ TM acting on a vector held V e j£(/ x E) gives 

TZ.V = TZ.(dt(V)dt + U,V) = dt(V)TZ.dt + TZ.Ii^V. (3.4) 

We can, hence, write TZ = Zdt + TZ with Z := TZ.dt and TZ := TZ.U^. Notice that 
TZ .V € X{M) is tangent to the submanifold Z('Et) for each tel. 

We can apply the preceding results regarding the decompositions of the metric 
and diheomorphisms to the particular case of g = Z*g. Let Vi, V 2 e X(/ x S) then 

{Z*g){Vi,V 2 ) = g{TZ.Vi,TZ.V 2 ) = g{dt{Vi)Z + T^.Vi,dt{V 2 )Z + TZ.V 2 ) 

= g{Z, Z)dt{V,)dt{V 2 ) + dtiV,)giZ,TZ.V 2 ) + dt{V 2 )g{Z,TZ.V,) 

+ g{TZyuTZy2). 

Comparing the previous expression with (3.2), we write 

7z = 9(2y-,ILZ-), (3.5) 

Pz = giZ,TZ-), (3.6) 

£iV| = g{Z,Z)-r\^zJz)- (3.7) 

Notice that 7 z is a symmetric (2,0)-tensor held on J x E but it is not a metric as 

dz{dt,-) = 0 . 

For each f € / let us dehne the map : E -> J x E : s 1 -^ (t, s), we have then that 
jt(E) = Tt and Zt := Z o : E M is an embedding for each tel. The fact that Z 
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is a diffeomorphism of the type considered in the paper guarantees that the images 
ZtCE), t ^ I, provide a foliation of M by spacelike hypersurfaces. The diffeomorphism 
itself can be reinterpreted as a smooth^ curve of embeddings. By using the jt map 
we can pullback the tensors (3.5)-(3.7) to S. This way we get 


IZt 

■= ftlz = Zpg. 

(3.8) 

hzuZt) 

:= ft^z = g{Zt,TZf), 

(3.9) 

eN'^ ■ 

^ {Zt,Zt) 

■- ^z - 9{Zt, Zt) - lz]ifd(Zt,Zt)^ f^(Zt,Zt)) ■ 

(3.10) 


Although in the next section we will discuss with more detail the manifold of em¬ 
beddings, it is important to highlight at this point that the Riemannian metric 7 ^^ 
depends only on Zt whereas both and ^ ^ depend also on the velocity 

Zt, which is a vector held in M along the embedding Zt- Moreover, is linear 

in Zt as a consequence of (3.9), as also is because 

^{ZuZt) = := eginz,, Zt ), (3.11) 

where we denote nzt the future-pointing unit ^f-normal to ^t(S) c M. In a similar 
fashion we introduce 9zt as the future-pointing unit normal to Zt{d'E) c dj^M. 

We introduce now two objects that will be useful in the following. If X : S ^ M 
is an embedding we dehne 


fe)^=(TX)^ (3.12) 

(ev)^=9„^(TX)f7J^ (3.13) 


Here and in the following we will use the Penrose abstract index notation when 
convenient (a, jS, etc. are abstract indices on M and a, b, etc. abstract indices on E). 
These objects satisfy 


(rx)a(ex)!; = <^a , (w)a(ex)^ = en'^rixp , 

where . From equation (3.9) we can write 

It is straightforward to write now (2.2) in the form 


(3.14) 

(3.15) 




/r.i /j.,. ^ M 5 -s''’ 


- fdT 

2 Jl JdT. 


\L 


'ipt-{d'ipt)b{eztfpZ^t 


nztiZt) 


eez,{Zt)bl ^vol 


IdZt 


+ d^Jt) - m^i^t I vol^^^ 

(3.16) 


where i>t ■= 'ipt ■= Jti’, 'ip ■= dii{dt), hz^ := BoZto iq, jgzt ■= ^llZt and : (9S S. 
Although it is possible to read off the formal expression of the Lagrangian for our 
model from the previous equation, it is important to pay attention to its domain, an 
appropriate subset of the tangent bundle of the conhguration space. We devote the 
next section to this issue. 


^We will not discuss topological issues here but just mention that they can be handled by the 
convenient calculus approach of Michor [23] . 
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4 Dynamics: geometric arena 

Configuration Space 

The configuration spaces of the different parametrized scalar fields systems that 
we consider in the following sections are 

Cd ■■= C'S’(S) X Emb^_,i(E,M) c x Emb,.,i(S,M), (4.1) 

Cr := X EmbJ,i(E,M) c C'-CS) x Emb,.,i(S,M). (4.2) 

Here (^^(S) = (^““(S) consists of smooth scalar fields while the elements of (7^(2) 
are smooth functions that vanish at the boundary 5E. Embg_si(E,M) is the space of 
smooth, spacelike embeddings of S in {M,g). This space has been studied in [3,4] 
(see also [24]). Our space Embg_si(S, M) is the subclass of such embeddings taking 5S 
to dj^M. For J 6 {D, R) we will denote a configuration {q,X) e Cj where g : S ^ M is 
a smooth scalar field and X : S ^ M is an embedding. 

Velocity Space 

In the following we will use the standard notation for fibre bundles so, for in¬ 
stance, r(X’"rM) will denote sections of the pullback bundle X*TM defined by the 
embedding X. The velocity phase-space TCj is the product manifold 

TCj = TCy{E) X rEmb^_3i(S, M), (4.3) 

where TC'J°(S) = Cy{T,) x Cy{T,) and for each X e Emhg_^i{T,, M) we have: 

TxEmbl,i{E,M) = T^(X*TM) := {Vx e r(X*TM) : Vx\dE e r(X’^TasM)}. 

This can be restated in terms of the commutativity of the following diagrams 

TM Td^M 



where the vertical arrows represent the natural projections in the respective tangent 
bundles. The condition defined by the second diagram implies 

g{Vx\dE,ux)=0, (4.4) 

where ux = V o X o iq and V is the unit, spacelike, outer normal to dj^M. For generic 
elements of rxEmbg_si(S, M) = r(X*rM) only the left diagram applies, as is the 
case for nx- In the following, a typical element of the velocity phase-space T(^q,x)Cj 
will be denoted as V(q,x) = {v, Vx), where v e Cy{E) and Vx e r^(X*TM). 

A comment about the notation that we have used so far and will be used through¬ 
out the paper is in order now: if we consider some tensor field V e r(r^’®Emb(S, M)) 






over the space of embeddings we will denote Vx e T^^Emb(S, M). The same subindex 
notation will be applied for more general objects like r and e, that whenever they 
are considered over a fixed embedding X, will be denoted txXx- Therefore if we 
consider V e X(Emb(S,M)) then Vx e TxEmbg_si(S, M), and we have just seen that 
this element can be considered as a vector held over the embedding X. It is clear 
then that such vector held can be decomposed in the form Vx = VJ^nx + tx-Vx^ where 
^Vx 9(yx,nx) e C'J’(S) and V^ e X(S) is dehned by tx-Vx ■= Vx - Vx'^x- Such 
decomposition can be made over the space of embeddings writing simply 

V = V^n + T.V^ e X(Emb(S, M)). 

Altough here both addens are vector helds on the space of embedding, notice that 
:Emb(S,M) ^C'°°(S) and EEEmb(S,M) ^ X(S). 

Finally notice that if M = / tangent bundle TEmbg_si(S, M) is the trivial 

bundle TEmbg_si(S, M) = Embg_si(S,M) x 


Phase Space 

The phase-space T*Cj is the product manifold 

T*Cj = T*Cy{E) X T*Emb3_,i(S,M), (4.5) 

where = Cj{T,) x C7(S)' and 

TJ(-Embg_si(S, M) = {Px \ Px • T^{X*TM) R linear and continuous} . (4.6) 

In the case M = / x R^ we have that T*Embg_si(S, M) is a trivial bundle over the 
base Embg_si(S, M). A typical point of is of the form P(q,x) = (P,Px) where 

p e C'J’(S)' is a distribution, i.e. a continuous linear functional p ■ C'J’(S) R, and 
Px ■ T^(X*TM) R is a continuous linear functional dehned on the space of vector 
helds along the embedding X. The phase-space T*Cj is equipped with the symplectic 
form n dehned by 

np<.,„(U,U) = U2(Ui) - Ui (y,2) + U.2(Ui) - Vpi(rx2), ( 4 . 7 ) 


where 


y := (((«, X), (p, Px)), ((Ui, U.). (Ui. UO)) e 7’p<.,x>r*0 ■ 

As Yi are tangent vectors of the phase space, we have Y^i e CJ’(S), lx* ^ r^(X*TM), 
Ypi e C'J°(S)' and Ypi : r'^(X*TM) R. Notice that the sans-serif subindices are 
used to denote the components of the vector held, in particular the subindex X of lx 
is entirely diherent from the subindex X used, for instance, in nx, where it denotes 
the evaluation of the vector held n € X(Emb(S,M)) at the point X. 

As we will see, the distributions p and Px that we will need in the paper can 
be dehned in terms of a scalar held p e (^““(S) and two covector helds along X; 
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Px ■■ T9{X*TM) ^ and Pax ■ T^iX^Tdj^M) -> C°°(aS). They have the form 


p(v) = / pnvols, (4.8) 

JT, 

PxiV)= [ Px(V)yoh+ [ PaxiV)vol9j:. (4.9) 

Js Jas 

Here vols and vol^s are hxed volnme forms on S and 9E, respectively. These distri- 
bntions can be alternatively expressed in terms of the metric volnme elements of the 
Riemannian metrics 7x = X*g and jax ■= {X o ia)*g = Pg'^x'- 


p{v) = / pnvols = / -^vol.^^ , 

Px{V)= f Px{V)voh+ f Pax(R)volas= f 
JT. JdT JT 





PaxiV) 


vol 


"fdx ■ 


After having introdnced these geometric elements we discnss now the Hamiltonian 
formnlation. 


5 Hamiltonian formulation 


From the 3 + 1 expression (3.16) of the action, we can dehne the Lagrangians 
Lj : ^ M for J € {D, R} given by 

^j(v(q,x)) = nx(R)^^ ) 

L ^^x{V)b^xQ‘^vo\^^x ■ ( 5 - 1 ) 

^ w/ 

where 


Pj ■= {v(q,x) £ TCj ■■ enx(V) > O} . 


The Hamiltonian formnlation of dynamics is dehned in the cotangent bnndle of 
the conhgnration space. In order to go from the tangent to the cotangent bnndle we 
need to compnte the hbre derivative dehned by onr Lagrangian. This is simply 


FLj{y(^q^x)) • V(q,x) - - 

1 



2 Jt 

-\L. 

'"It 


n-ex(dg,y) 

enxiy) 

e'y]l{dq,dq) - rri^q^ - 


ex(dg,R) vol^^ 


' v-ex{dqy) ' 
. nxiV) , 


(5.2) 


\enx{V) vol 


lx 


e9x{y)hxq^xol 


IdX 


«-ex(iq,V) ^ 

enxiV) 
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The fibre derivative defines, among other things, the canonical momenta. Calling 


{Px)a - 

s/lx 


v-ex{<^qy)' 


{exTMq)a 


enx{V) 
e'y/-(dq,dq) - rr/q^ - 


' v-ex{dqy)' 

, nxiy) 


£{nx)a , 


yid. 


{Pax)a = -^bWe{ex)a, 

v-eAd^ — 

^ en^(V) 


it can be seen that the distribntions defined in (5.2) belong to the class defined by 
eqnations (4.8)-(4.9). If we introdnce now 


Pa{q,p) = p{dq)a, (5.3) 

P-i{q,P,X) = ^^^{ryq‘^ + — -e-f]l{dq,dq)], (5.4) 

2 V 7x / 

K(g,X) = lybW\8^. ( 5 - 5 ) 

the primary constraint snbmanifold Jv['l := FLj(T>j) can be written as Adf = {Ca = 
Oa,Ca = ^a} wliere 


cyq,x,p,p) := {Px)a +£inx)a'Hi{q,P,l() + {ex)aPa{q,p), (5.6) 

C^{q,X,Pa) ■■= {Pax)a + e{ex)anf{q,X). (5.7) 


As the Lagrangian is homogeneons of degree one, it is immediate that the energy 

E{^ig,X)) = FLj{^{q,X)) ■ ^(q,X) “ ^j(V(g,X)) (5.8) 

vanishes in the fnll tangent bundle of the configuration space. The Hamiltonian H is 
defined on Aif as E =: H o FLj and, hence, it is also zero. 

The equation defining the Hamiltonian vector fields Z on A4f is simply 


izu; = 0 (5.9) 

where u is the pullback of H to M.( and Z := (({q,X), (p, P)), ((.^q, Zx), (Zp, Zp))) 
denotes a vector field on A4f. Notice that the components (Zp,Zp) are of the form 
(4.8) and (4.9) respectively. 

The easiest way to solve the equation for the Hamiltonian vector field is to find 
the solutions to u)(Y,Z) = 0 for all Y e As we are looking at as a 

submanifold embedded in phase space we must find a convenient way to characterize 
and use these fields. In the present case this can be done by requiring the vanishing 
of the directional derivatives of the constraints along the field. 
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Let us consider a point V{q,x) - (Py PyPd)) in phase space, then oj is given 

by 

‘"p,„x,(VZ) = ^ - nz, + ZppU" - rp<.Z?)vol,, ( 5 . 10 ) 

JdT. V7ax 

where 

Yp„ = L> ^ - eua'Hi - , Zpa = d{^- euaUi-y ( 5 - 11 ) 

Y^^ = D[-eejii) , Z^^ = D[-eeaHf) (5.12) 

denote the variations (A.l) and (A.3) given in the appendix A. The concrete form of 
u(Y,Z), obtained by a long but straightforward computation, is given by 


u{Y,Z) = 


L 

Z,-( 

-L 


L 

n-( 

L 

Z,-( 


-^en^ + e^^{dq)bj Z^ 


_p_ 

v/7 


ena +ei{dq)b]Y^ 



Wx^grad.^g - j 

/ \ ~ 

—^ + nlqZy + div^Y 
Lx /7 "" 

^eZ^gradp - j 


(5.13) 

VOl.y 

voL 


+ el^{dq)bjY^ 
^£na + e^((ig)bjzx 


eZ. 


1 / P 


,x/7 


en/s + P^fj{dq)a - eb^qu/s ppjvola.^ 


u' 


sY^ en /3 + e‘p{dq)a - eb'^qu/sj j^vola.^ , 


where we have written ux = 'r-zz]^ + iz-‘-rix, Wx\^ •= 7v(i^xi have dropped the 

X subindex in all the elements to render the last expression more compact. Notice 
that due to the antisymmetry of w, the Weingarten maps Kx and appearing in 
equations (A.l) and (A.3), cancel out in the hnal form of (5.13). 

Requiring the previous equation to be zero for every Y e we hnd that the 

Hamiltonian vector held in the interior of S is 




— ;- Zv + Cz^Q ~ 

VTx 


P 


,\/W 


enxa + {ex)a{dq)o 


-rri^V^ gZx + div 


P 


lx 


.'/ir 


ZJ-sZ,;grad ^^9 


(5.14) 

(5.15) 


for all the types of boundary conditions that we consider here and, also, when S has 
no boundary. Notice that, as we are working with smooth objects, equations (5.14, 
5.15) must be extended by continuity to 9S. In the preceding equations the gradient 
and divergence are dehned in terms of the metric 7 ^ in the standard way, Cy^ denotes 
the Lie derivative along a vector held € X(S) and we have used the decomposition 


Zx = Z^n + r.Zx . 
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When solving equation (5.9) no restriction over Zx in the interior of S arises, hence 
it can be chosen freely (within the class of regular objects that we are considering). 

Equations (5.14, 5.15) coincide with those appearing in the classic papers on the 
subject [3]. For closed spatial manifolds this result is all that is needed to get the full 
Hamiltonian description, however, in the presence of boundaries some extra conditions 
may appear. In order to see this notice that, although the last boundary integral in 
(5.13) vanishes, because (5.14) is also true at the boundary, we need to require that 



P 


—en„ + e”(dg)fejyx" 




P 


—en/s + e«(dg)a - eb^qu/s 
V7 




\iy^\ ’ 


(5.16) 


be zero for every Y € 


5.1 Dirichlet boundary conditions 

As we mentioned at the beginning of the paper, the way to deal with homogeneous 
Dirichlet and Robin boundary conditions differs in some important details. Dirichlet 
boundary conditions are enforced by restricting the configuration space to those scalar 
fields which vanish at the boundary. As a consequence of this, the q component 
of the admissible vector fields Y e j£(A4f) must also vanish at the boundary, in 
particular the q component (5.14) of the Hamiltonian vector field must be zero. By 
continuity this leads to 
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Zx, 


——enx + Tx.grad g 

Vtx 


as 


—=enxa + (ex)Ji(dg)a 


as = 0. 


(5.17) 


We can check now that this expression implies that the boundary integral term (5.16) 
vanishes and, hence, there are no extra conditions on the vector held. Indeed, hrst 
notice that as .^x|as e r^(X*Td^M), whenever is not zero at the boundary aS, 
equation (5.17) implies 


P 


enx + Tx-grad^^q oc ux , 


(5.18) 


and then, the factor 




—=enxa + (ex)i{dq)b 

VTx 


\za 


(5.19) 


of (5.16) is zero, because Hq|aS = 0 and g(Yx,u)\dT, = 0 as a consequence of (4.4). 

It is useful to rewrite equation (5.18) in a different way. If we denote 2;x •= -^xias = 
z^n + T.z^, the condition Zx\d'Z € T^{X*TdY.M) becomes 


£z/j^4 + 7x(z/]f,4) = 0, 


(5.20) 


while equation (5.17) is equivalent to 


P 1 
\/lx 


+ 7x(grad.^^g,4) = 0 


at as. 


(5.21) 
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Equations (5.20) and (5.21) define an homogeneous linear system for the boundary 
values that has to be solved in terms of p and dg. We always have ^ 

Zx = 0 as solutions. Pointwise, a non-trivial solution is only possible if 


eiz^grad q 


P 




-u 


X 


as 


(5.22) 


is zero at s e aS. These conditions are not constraints in the standard sense but rather 
dehne sectors in the primary constraint manifold. The best setting to understand 
the appearance of these sectors is the study of the parametrized electromagnetic 
held (see [20]) where the all-important Gauss law appears precisely in the same way. 
Labelling the points of by the support of (5.22), we characterize the relevant 
sectors where the “number and type” of the independent components of zx changes. A 
dinamically relevant sector is the one with elements {q,p,X) such that the associated 
expression given by (5.22) has empty support (or equivalently the condition (5.18) 
holds everywhere). In that case zx can be different from zero, and so the dynamics of 
the parametrized system allows the embeddings “to advance in time” and generate a 
genuine (local) spacetime foliation. If, instead, the zx were zero the integral curve of 
embeddings would be forced to be “stuck” to the same section Ao(S) of d^M and, 
hence, would not be suitable to describe held dynamics in the usual way. 

A consistency condition must be imposed now: the Hamiltonian evolution dehned 
by the Hamiltonian vector helds that we have obtained must be compatible with 
condition (5.18) (equivalently the vanishing of (5.22)), which is a requirement to 
have proper dynamics for the embeddings as we have just seen. This means that 
the Hamiltonian vector helds must be tangent to the submanifold formed by the 
elements {q,p, X) for which (5.18) holds. This points out to the existence of additional 
requirements necessary to have consistent dynamics. This situation exactly mimics 
the one found in the Hamiltonian treatment of the scalar held with Dirichlet boundary 
conditions (see [21]) and is in perfect agreement with the known results for the scalar 
held in the smooth case [22]. From here on, the determination of the inhnite chain of 
conditions necessary to have well dehned dynamics for smooth held and embeddings, 
follows exactly the steps of the GNH algorithm as the main geometric issue involved 
is the tangency of the Hamiltonian vector helds to the submanifolds dehned by the 
successive conditions. Although they can be obtained in a straightforward manner 
using Table 1 of appendix A, they are somewhat complicated and their particular 
form is not specially illuminating, hence we do not give them here. 

As a hnal remark notice that if = 0 at X(dT,), or equivalently is tangent 
to X(E) at A(5E), the condition (5.22) immediately implies p|9E = 0 (as never 
vanishes). This condition is found in the Hamiltonian treatment of the scalar held 
with Dirichlet boundary conditions [21], where a special foliation with Zx = 1 and 
Zx = 0 is used. 


5.2 Robin boundary conditions 

In order to implement the Robin boundary conditions we allow the helds q to 
take non-zero values at 9E and consider b arbitrary (if it is zero we would be in the 
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Neumann case). In general the boundary integral (5.16) would not vanish as before, 
so according to equation (5.9) we have to require it to be zero for every Y e X(A^f). 
Denoting again zx == we have 


P 


_s/lx 


enx0 + {exYfi{M)a - £h\qvxp 


'X 


as = 0 


(5.23) 


which plays a role analogous to the one of equation (5.17) in the Dirichlet case, despite 
their very different origin. 

Wherever t 0, equation (5.23) implies 


P 

\/7x 


enx + Tx-grad^^q - eh\qvx e Txat^dj^M. 


(5.24) 


Condition (5.24) is the analogous of (5.18). It is interesting to mention that in the 
Dirichlet case we have that 


——enx + Tx.grad g 

s/lY 

is normal to d^M, in the Neumann case {b = 0) it is tangent to it, and in the Robin 
it is neither normal nor tangent. Clearly we can recover the Neumann case from the 
Robin one by making 6 = 0, however, there is no such way to pass from Robin to 
Dirichlet. This shows once more the intrinsically different nature of these boundary 
conditions. 


As in the Dirichlet case, we can reexpress (5.24) in a language based entirely on 
objects dehned on S. We have the equations 


P 




+ dq{n\) - eh\q 


as = o 


+ = 0 - 


(5.25) 

(5.26) 


Again, these should be considered as a system of equations for which has 

non-trivial solutions when the following compatibility requirement holds 

on cIS. 


2 „ P 


dq{n\) = ebxq- —=>yx 

VTx 


(5.27) 


Notice that the vector held 






YixWxyx) 

is the unit outer normal (according to yx) to and, hence, 

dg(i^i) = Yixyx^^x) 

is proportional to the normal derivative of the scalar held q. 
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We can see again the same phenomenon that we fonnd for the Dirichlet bonndary 
conditions. There are sectors in the primary constraint hypersnrface Adf where the 
components zy^ mnst necessarily vanish whereas in other parts of Adf, formed by 
elements {q,p,X) for which (5.24)/(5.27) hold, the zy components may be different 
from zero. 

The sitnation at this point is conceptnally eqnivalent to the one that we explained 
in the Dirichlet case: the conditions (5.24)/(5.27) behave as secondary constraints 
whose stability nnder the dynamics given by the Hamiltonian vector helds mnst be 
enforced. From a geometric point of view this is a tangency reqnirement that pro¬ 
vides an inhnite chain of constraints (whose relatives in the nnparametrized case are 
explicitly given in [21]). 

We can make the same hnal remark as in the Dirichlet case, if we consider embed¬ 
dings X snch that the normal to the bonndary ux and the normal nx to ^(S) are 
orthogonal (i.e. snch that =0), and taking into acconnt that, in this case, = Cx, 
we recover the condition 

[dqi^x)-ebj,q)\dJ: = 0 (5.28) 

that we fonnd for the non-parametrized case in [21]. 

6 Comments and conclusions 

The main resnlt of the paper is the precise description of the Hamiltonian formula¬ 
tion for a parametrized scalar held dehned in a bounded spatial region with or without 
boundaries. When no boundaries are present we recover the results of [3]. We have 
considered boundary conditions of the Dirichlet and Robin types (Neumann bound¬ 
ary conditions are a particular example of the latter) and worked within the class of 
smooth helds and spacelike embeddings. We have obtained the constraint submani¬ 
fold in phase space and the concrete form of the Hamiltonian vector helds. Contrary 
to the statements in [8], we have found no obstructions to get the full description 
in the Robin case. We expect that our detailed formulation can be used to extend 
the polymer quantization of the scalar held to spatial manifold with boundaries, in 
particular in the 1-1-1 dimensional case. Other systems that can be considered from 
this perspective are Einstein-Rosen waves coupled with massless scalar helds [10] 

There are a number of interesting facts that we would like to mention. First, the 
solutions to equation izOJ = 0, giving the Hamiltonian vector helds, can be obtained 
pointwise on the primary constraint hypersnrface, however, it is a non-trivial issue 
to understand to what extent they dehne smooth vector helds on Xi'l or a submani¬ 
fold thereof. There are a number of regularity issues that must be characterized and 
understood. For instance, there are sectors in the primary constraint hypersnrface 
where some components of the vector helds are forced to be zero whereas in others, 
dehned by the vanishing of specihc functions of the conhguration variables and mo¬ 
menta, they can be diherent from zero. We have paid particular attention to a sector 
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which is specially important from the dynamical point of view: the one dehned by 
conditions (5.18) and (5.24) for the Dirichlet and Robin cases respectively. 

Second, for parametrized models the energy is zero and, hence, the appearance of 
different sectors is even more important in the Lagrangian symplectic approach. In 
this approach the canonical symplectic structure is pulled back to the tangent bundle 
of the conhguration space with the help of the hbre derivative, and the dynamics is 
obtained by hnding the Hamiltonian vector helds given by its degenerate directions 
[16]. The standard Hamiltonian approach can be mirrored in this setting, where no 
constraints show up and all the subtleties associated with the singularity of the system 
must manifest themselves under the guise of regularity issues of a nature similar to 
the one of those discussed in this paper. 

Finally, we would like to point out that the Hamiltonian description accommo¬ 
dates, in a straightforward way, some non stationary situations where the metric qq 
induced on the boundary d-^M has no timelike Killing vector helds. In the Robin case, 
it is furthermore possible to introduce non-constant boundary conditions encoded on 
the value of the scalar function B on dj^M. These non-stationary situations ht quite 
naturally in our approach as stationarity plays no role in any part of it but, of course, 
the integrability of the Hamiltonian vector helds must be considered with due care. 

The method that we have used in the paper is quite general and can be applied in 
a similar way to other parametrized held theories (gauge helds, gravity, etc.) dehned 
in diherent types of spatial regions (bounded and/or unbounded). We plan to explore 
them in the near future. 
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Appendix A Variations in Emb(I],A^) 

A detailed description of the inhnite dimensional manifold Emb(S,M) of embed¬ 
dings X : S ^ M in the case where {M,g) is a Riemannian manifold can be found 
in [24]. References [3,4,12] provide a summary on the relevant results concerning the 
manifold Embg_si(S, M). We refer the reader to those papers and references therein 
for further details. In this appendix we just list a number of expressions for the 
variations of several geometric tensor helds on Emb(S,M) that are necessary to get 
the Hamiltonian formulation for the systems that we consider here. We give them in 
Table 1, where we have corrected the typos in equations (9.8) and (9.10) appearing 
in section 9 of [12]. 
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fx 

^(X,X)f 

bx-=BoX 

7x-,x)'> = dB(X) 

vx-=V°X 

D(xX)^ = Yx^ 

(ev)S 

= (KtX^^ + ej^\dX^),) en^ - - X^K^ 


D(X.X,^° = pfai,V™ + £(d^0«) en“ + T,“ (V.X" - 


C(x,x)"“ = -C(a7A-’‘ + 

(^A')a ■” QajS'^x 

= -e“ (XabX"' + £(dX")„) 

i'yx)ab 

“ dfecVaV^'^ + 7acVfeX^'^ - 2X^Kab 

7x 

+ 2X^K'^’’ 

vol^^ 

^(X,X)V0l-y = (Va^"“ - X^K:)vo\^ 


7x.x)'/7 = hT~-vvc)v^ 


Table 1: Variations of the relevant geometric tensor fields on Emb(S,M). Here and 
in the following, denotes the Weingarten map associated with V(E) c M, with 
indices lowered and raised with the help of jab and 7 “^, and V denotes the Levi-Civita 
connection of (^, 7 ). 


Using Table 1 , it is immediate to obtain Ypa = D {-enaV-i - in the form 

(dropping the X snbindex) 


v/7 


= - (+ e|((dg);,) Yp + 1 7 “‘=(dg)cen„ - e“) (dyq)„ - erio^m^qY^ 


,v/7 


V7; 


-^(dg), + er\^Y^)a) - e (y^YX’ - (dg)„(dg), 




,\/7 7 


eria 


(AT) 


+ e“ 


n, 


+ <!(drx‘)4 


Snbstitnting the ambient spaces M and S by and and considering the 
indnced embeddings dj:,M we can obtain a similar table of variations, where 

now the relevant objects are associated with the bonndary 9S. For instance we have 
to nse the Weingarten map of X{d'E) c d-^M and the Levi-Civita connection 
of {dTj,'yg). Also, for any given Wx e r^(A*TM), we have to nse the decomposition 


Wx\dE = W^Ox + 


(A. 2 ) 


18 















where := TX®, := g{Wx\dT,,9x) and e X(9S) is defined by ^x-^x 

Wx-W^^9x- With all these elements we find that Yp^ = D {-e9aT-Lf) can be expressed 
as (dropping again the X subindex) 


■/Ts 


^^aydja _ + h (Y^^9^{dh)p + Y^'^P{dh)p) + h'^qY^ 


Xf,Ff^ + £(dyx"), 


e9^ 


(A.3) 


As we have to make contact with the decomposition W = W^n + t.W^ associated 
with n and r, we need to use the following relations 




Vpi 


= ppT ^ 






T|2 


-V 


(A.4) 


Plugging (A.4) in (A.3) and using (A.l), we obtain the final form of (5.13). 

An interesting application of the results of Table 1 is the derivation of the hyper¬ 
surface deformation algebra [12]. Considering V = C-'-n + r.W and W = IP-'-n + r.VW, 
vector fields on the space of embeddings, we can define 


[C,W] :=DvW-DwV 


or, equivalently, [C, hP]x := D(^x,Vx)^ ~ D(^x,Wx)^■ Froni these, it is immediate to 
obtain the decomposition 

[C, VP] = (DyVP^ - DwV^ + dV^{W^) - dVP^(W))n 

+ r.(DyVP" - DwV^ + e(Wgrad.^^VP" - VP^grad.^^ W) - [W, VP"]) . 
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